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Null Energy Condition (NEC) requires the equation of state (EoS) of the universe w u satisfy 
w u > — 1, which implies, for instance in a universe with matter and dark energy dominating w u — 
Wm^lm + Wde^ide = w decide > — 1. In this paper we study constraints on the dark energy models 
from the requirement of the NEC. We will show that with ilde ~ 0.7, Wde < —1 at present epoch is 
possible. However, NEC excludes the possibility of uide < — 1 forever as happened in the Phantom 
model, but if Wd e < —1 stays for a short period of time as predicted in the Quintom theory 
NEC can be satisfied. We take three examples of Quintom models of dark energy, namely the 
phenomenological EoS, the two-scalar-field model and the single scalar model with a modified Dirac- 
Born-Infeld (DBI) lagrangian to show how this happens. 
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It is well known that energy conditions play an important role in classical theory of general relativity and 
thermodynamics In classical general relativity it is usually convenient and efficient to restrict a physical sys- 
tem to satisfy one or some of energy conditions for study, for example, in the proof of Hawking-Penrose singularity 
theorem^ [U, the positive mass theorem^] and so on, while in thermodynamics energy conditions are the bases for 
obtaining entropy bounds p, [y]. Among those energy conditions, the null energy condition is the weakest one which 
states that for any null vector n^ 1 the stress energy tensor T av should satisfy the relation 
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In general, the violation of NEC leads to the breakdown of causality in general relativity and the violation of 
the second law of thermodynamics [7J. These pathologies require that the total stress tensor in a physical spacetime 
manifold should obey the NEC. In the framework of the standard 4-dimensional Friedmann- Robertson- Walker (FRW) 
cosmology the NEC implies p + p > 0, which in turn gives rise to a constraint on the equation of state of the universe 
(EoS) w u defined as the ratio of pressure to energy density, w u > — 1. In this paper we study the constraints on the 
dark energy models from the requirement of w u > — 1. 

In the early Universe with radiation dominant the EoS of the universe w u is approximately equal to ^ and in the 
matter dominant period w u is nearly zero, so NEC is satisfied well. However when the dark energy component is not 
negligible we have 
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w u = w m VL m + w de fl de > -f , (2) 
where the subscripts 'm' and 'cfe' stand for matter and dark energy, respectively. With w m = 0, inequality @ becomes 

Wde^de > -1 • (3) 
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From the inequality above, we can see that models of dark energy with Wde > — 1 such as the Cosmological Constant 
and the Quintessence satisfy the NEC, but the models with Wde < — 1 predicted for instance by the Phantom theory 
where the kinetic term of the scalar field has a wrong sign does not. Interestingly we can see that NEC might be 
satisfied in models if w de < — I stays for a short period of time during the evolution of the universe. In this paper we 
will show this happens in the Quintom models of dark energy. 

The paper is organized as follows: in section II we will present three examples of the Quintom models to show how 
the NEC is satisfied and the section III is the summary of the paper. 



II. NULL ENERGY CONDITION AND THE QUINTOM DARK ENERGY 

Quintom is a dynamical model of dark energy Q. It differs from the Cosmological Constant, Quintessence @, 
Phantom K-essence[ll| and so on in the determination of the cosmological evolution. Although the current data 
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in combination with the 3-year WMAP[l2|, the recently released 182 SNIa Gold sampled! and also other cosmological 
observational data show the consistence of the Cosmological Constant, it is worth noting that the dynamical dark 
energy models are not excluded and Quintom dark energy is mildly favored (for recent references see e.g. 
The most salient feature of the Quintom model is that its EoS can smoothly cross — 1. In this section we will study 
the implications of NEC on the Quintom models. Working with three specific examples we will show the NEC can 
indeed be satisfied. 



A. A phenomenological model with parameterized EoS across — 1 

With a simple calculation, the inequality ([3]) can be rewritten as 

(l + w de (a))n de0 exp^ [-3(l + ifl« fa (a / ))]dWj+^>0 1 (4) 

where the subscript '0' represents today's value. 

We firstly start with a phenomenological model with a parameterized EoS which will be able to cross over — 1 : 

w d e = w +wi(l - a) , (5) 

where a is the scale factor which we normalize to be ao = 1 at present. One can see that when a is equal to 
(1 + wo + wi)/wi, the EoS of dark energy crosses — 1. This type of parametrization for the EoS has been widely used 
in the literature [T3, for the fitting of constraining the EoS to the observational data. 

Inserting Eq. into the inequality (0]), the NEC puts a constraint on the parameters wq and w\. In Fig. [TJ we 
take wq = —0.9 and w\ — —0.3, and then plot the evolution of the dark energy EoS and also the EoS of the universe 
respectively. One can see from Fig. [T]that the EoS of dark energy Wde crosses —1 from below to above. In this case 
the EoS of the universe w u evolves from zero which corresponds to the matter dominant epoch, and then reaches its 
minimal value during which the universe enters the dark energy dominant period, and finally returns to zero in the 
future. We can read directly from this figure that the minimal value of w u stays above —1 which satisfies the NEC. 
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Figure 1: Plot of the evolution of the EoS Wde with Wde = wo + wi(l — a), and the EoS w u for the universe as a function of 
In a. Here in the numerical calculation we have taken wo — —0.9 and mi = —0.3. 

One could take another example with a different parametrization of the EoS: 
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w de = wo + wx(l - a + —a 2 ) . (6) 
lo 
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This example is different from the previous one since we introduce a square term of the scale factor which makes the 
EoS Wde crosses —1 twice. Taking proper values of those parameters, this model can give a scenario that dark energy 
stays in the Phantom-like state only for a while and then returns to be Quintessence-like. After taking the similar 
calculation, we plot the evolution of the EoS of the model and the universe respectively in Fig. O In this example the 
EoS of dark energy Wde starts to evolve from Quintessence-like in the past, and then enters the Phantom- like state 
for a short period of time, and finally returns to above —1. We also find that w u behaves similar to that in Fig. [TJ Its 
value resides on zero in the past and in the future which means that matter have dominated and will dominate again 
the evolution of the universe. Although w u runs away from zero during the evolution, its value keeps being larger 
than —1 as shown in Fig. [2l Therefore, this model can be consistent with the NEC as well. 
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Figure 2: Plot of the evolution of the EoS Wde with Wde = Wo +Wi (1 — a + j^a 2 ), and the EoS w u for the universe as a function 
of In a. Here in the numerical calculation we have taken wo = —0.93 and wi = 0.23. 



B. The Two-scalar- field Quintom model 



Building the field model ofQuintom dark energy is a challenge due to the No-Go theorem which has been proved 
in Ref. [l9|(also see Ref. @, [13, ISI, [H, H [H). This No-Go theorem forbids a traditional scalar field model with a 
lagrangian of general form C — £(<fi, V M (/>VV) from having its EoS cross over the cosmological constant boundary. 
According to this theorem, dynamical models like Quintessence, Phantom and K-essence are unable to realize their 
EoS cross — 1. Therefore, to realize a viable Quintom field model in the framework of Einstein's gravity theory, 
one needs to introduce extra degrees of freedom to the conventional theory with a single scalar field. The simplest 
Quintom model is constructed by two scalars with one being Quintessence-like and another Phantom-like proposed 
firstly in Ref. [|, and this model has been widely studied in detail later on. In recent years there have been a lot of 
activities in the theoretical study on building Quintom models, such as a single scalar with high-derivative l2ol |26|. 
vector field [27j. extended theory of gravity [28j. Lorentz-violating dark energy models [29j and so on, see e.g. [3d| |. 

In this section we investigate a two-scalar-field Quintom model of dark energy in flat FRW cosmology which is 
described by the action 



S = 



R 
16ttG 



de 



(7) 
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with 



where R is the Ricci scalar of the universe, C m is the lagrangian of matter, Cde is the lagrangian of dark energy, and 
the metric is in form of (+, — , — , — ). Here the field 4>i has a canonical kinetic term, but (f>2 is a ghost field. With Cde 
in |[8]), we can easily obtain the energy density and the pressure of this model, 

Pde = \$. -\&+V, (9) 

Pde = \i>\ - \<& - V , (10) 

where the dot denotes the derivative with respect to the cosmic time, and by the variational principle the Einstein 
equations are given by 

n &TrG . 1 • o 1 • o , , 

H 2 = —(^i-^i + v + Pm ), (ii) 

01+3^ + ^1=0, (12) 

2 + 3i70 2 -^1=0. (13) 

Since it is required that the total EoS of the universe satisfies the NEC, it is clear that the Phantom component of 
this kind of model can only be dominant for a short while and hence it gives a constraint on choosing the potential of 
the Phantom field. In the following numerical calculations, we will choose an appropriate potential of Phantom field 
and different potentials of Quintessence field, which can satisfy the NEC and give various Quintom scenarios. 

In Fig. [31 we choose the potentials of Quintessence and Phantom field to be Vfa = V\ exp(Ai</>f /M 2 ) and V^ 2 = 
V% exp(A202/^ 2 ): respectively. Here we normalize the fields with Planck scale and choose M to be O.OlAfp/. 1 During 
the evolution, the Phantom field will climb up along its potential and become dominant because its energy density 
will increase while that of Quintessence field will decrease. However, since there is a maximum of the potential, after 
reaching that point, the Phantom field will stay on the top and stop affecting the evolvement of the universe. From 
the figure, we can see that after the EoS of dark energy being less than —1 today for a short while, it can exit in the 
future and actually approach to —1. In this case we read from the figure that the EoS of the universe can be always 
larger than —1, and thus the NEC can be easily satisfied. 

In Fig. [5J we choose another potential form of the Quintessence field to be Vfa = hrn?<f)\ with the Phantom's 
potential unchanged. We can see from the figure that, similar to the reason above, Phantom can only dominate 
the universe for a while, and due to different potential of Quintessence, the EoS of dark energy model behaves 
very differently. The figure shows that after the universe exit the Phantom dominating phase, the Quintessence will 
dominate the universe again, and the EoS of dark energy will cross — 1 twice. In the future, however, this EoS will also 
become de-sitter like. In the whole process, the EoS of the universe will always be above the cosmological constant 
boundary, and this model of Quintom dark energy also satisfies the NEC. 

C. A single scalar field with a modified DBI lagrangian 

Having presented the examples of two-scalar- field Quintom models in consistent with the NEC, we in this section 
consider a class of Quintom models described by an effective lagrangian involving higher derivative operators. 

Due to the contribution of higher derivative terms, this kind of models can give rise to an EoS across —1 as pointed 
out in Ref . . A connection of this type of Quintom theory to the string theory has been considered in Ref . [3l| and 
. Here we take the string-inspired model in [3l[ for a detailed study to check whether it satisfies the requirement 
of the NEC. The action of this Quintom dark energy is given by 



'tie 



V{</>)y/l - a'V M 0V^ + l3'cj)U<j) . (14) 



1 Here and in the following, we would like to redefine the parameter Ai and A2 to be A1/M 2 and A2/M 2 
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Figure 3: Plot of the evolution of the EoS Wde for the two-scalar-field Quintom dark energy with the potential V t j >1 = 
Viexp(Ai0 2 ), V0 2 = V2 exp(A20f)i an d the EoS w u for the universe as a function of In a. Here in the numerical calcula- 
tion we have taken Vi = (3.1623eV) 4 , V 2 = (1.8612eV) 4 , Ai = -0.8 x 10 4 , A 2 = -0.5 x 10 4 and the initial values are fai = 0.02, 
<j>u = 5.0 x 10" 62 , cf> 2i = -0.01 and 4> 2i = 2.0 x 10^ 62 . 




Figure 4: Plot of the evolution of the EoS u>d e for the two-scalar-field Quintom dark energy with the potential V < / >1 = 
^m 2 (j)i, Vfe = V2 exp(A202)i an d the EoS w u for the universe as a function of lna. Here in the numerical calculation we 
have taken m = 1.0 x 10" 29 eV, V 2 = (1.8612eV) 4 , A 2 = -0.5 x 10 4 and the initial values are <j> u = 0.02, j> u = 5.0 x 10" 62 , 
<f> 2i = -0.01 and <j> 2i = 2.0 x 10" 62 . 
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This is a generalized version of "Born-Infeld" action[32|, [33| with the introduction of the 8' term(see, for example, 
[34| for derivation of the higher derivative term from string theory). To the lowest order, the Box-operator term cf>n6 
is equivalent to the term V^V^ when the tachyon is on the top of its potential. However when the tachyon rolls 
down from the top of the potential, these two terms exhibit different dynamical behavior. The two parameters a 1 and 



8' in (fT4|) could be arbitrary in the case of the background flux being turned on [35| . One interesting feature of this 
model is that it provides the possibility of its EoS Wde running across the cosmological constant boundary. In the 
following studies to make two parameters (a', 8') dimensionless, it is convenient to redefine a = a'M 4 and 8 = 8'M 4 
where M is an energy scale of the effective theory of the field. 
From (fl4|) we obtain the equation of motion for the scalar field 0: 

| D( ^ )+a V^, + M^/ + f^ = 0, (15, 

where / = — a'V ^4>V^6 + 8'60<fi and — dV/d<fi. Correspondingly, the energy stress tensor of Quintom dark 
energy is given by 



2M 4 py f ^ n M 4 f ^ UT 2M 4 ^ f 
2M 4 A f 



+^V„(^-)V M . (16) 



In order to simplify the calculation, we technically define another parameter ip = = — Jjjpj to solve (fT5j) and 
(fTB)) . In the framework of a flat FRW universe filled with a homogeneous scalar field 6, we have the equations of 
motion in forms of 

W M 4 V> V 2 . 86 2a ■ ■ 



and the energy density and the pressure of this field can be written as 

Pde = -^ -^-4M^ V - -ft' (19) 
onb -n ■■ 06 o M 4 ib 

According to the restriction of the NEC, we need this Quintom model to satisfy the inequality ([4|. Although in 
this case the phase space of Quintom dark energy is constrained, we will show below that the inequality Q can be 
satisfied easily. In the numerical study, we constrain the parameters a and 8 so that when expanding the derivative 
terms in the square root to the lowest order the model in (fll| gives rise to a canonical kinetic term for the scalar field 
6 [H|, i.e., a + P>0. 

Through calculating Eqs. (|17[) . (|18[) and the Fricdmann equations numerically and then comparing the results in 
the inequality ((4|), we can judge whether a model is consistent with the NEC. We first consider a model with the 
potential V{6) — cxp (_^^j. cxp (^0) ■ I n the calculation we normalize the fields 6 and iji with the Planck scale and 
choose the energy scale M — 1.2211 x 10 _4 ey. For a given set of the model parameters, we make the numerical 
calculations and plot the evolutions of the EoS of the Quintom dark energy and the universe, which are shown in Fig. 
[5j We can read that the EoS of dark energy Wde starts from a value larger than —1 then evolves to less than — 1, 
and soon exits the Phantom-like state, and eventually approaches the cosmological constant asymptotically. During 
the evolution the period of time for dark energy to be Phantom-like is very short. Therefore, the corresponding EoS 
of the universe evolves from the matter dominant period with w u — to the dark energy dominant period smoothly 
without violating the NEC as shown in Fig. [5] 

For another example we take a different form of potential V(6) = ^m?6 2 . Similar to what have been done above, 
we obtain another example consistent with the NEC. The normalization in this case is the same as that in the above 
one. In Fig. [51 we give the evolutions of the EoS of this Quintom model and the universe respectively. One can see 
from this figure the w u satisfies the NEC, but the detailed evolution of the universe differs from the one shown in 
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Figure 5: Plot of the evolution of the EoS Wde for the string-inspired Quintom dark energy with the potential V((j>) = 
. — r-X®. T-r-n- and the EoS w u for the universe as a function of In a. Here in the numerical calculation we have taken 

exp( — A0) + cxp(A0) 

A = 1CT S , V = (1.3183 x l(T 3 eV) 4 , a = 0.6, = 0.4 and the initial values are & = 12, $ = -2.3 x 10~ 5 , Vi = -3.84 x 10 2 
and ipi = —40. In choosing the initial values the prime represents the derivative with respect to In a. 



Fig. O In Fig. ©the EoS of Quintom dark energy Wde crosses the cosmological constant boundary from a fixed value 
below —1 to above and then evolves close to —1 in the future. In this case the EoS of the universe w u runs from zero 
when the universe is dominated by the matter, and then approaches the de-Sitter phase asymptotically. 



III. CONCLUSION 



In this paper we have studied the implications of NEC in the models of dark energy. We show that NEC excludes 
the models with Wde < — 1 forever as predicted by the Phantom dark energy, however allows the possibility of having 
Wde < — 1 for a short period of time as it happens in the Quintom models. We have shown explicitly in this paper 
three examples of Quintom models where NEC is satisfied. 
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